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Abstract. From the exact single step evolution equation of the two-point correlation 
function of a particle distribution subjected to a stochastic displacement field u(x), we 
derive different dynamical regimes when u(x) is iterated to build a velocity field. First 
we show that spatially uncorrelated fields u(x) lead to both standard and anomalous 
diffusion equation. When the field u(x) is spatially correlated each particle performs 
a simple free Brownian motion, but the trajectories of different particles result to be 
mutually correlated. The two-point statistical properties of the field u(x) induce two- 
point spatial correlations in the particle distribution satisfying a simple but non-trivial 
diffusion-like equation. These displacement-displacement correlations lead the system 
to three possible regimes: coalescence, simple clustering and a combination of the 
two. The existence of these different regimes, in the one-dimensional system, is shown 
through computer simulations and a simple theoretical argument. 
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1. Introduction 

The study of diffusion phenomena has a long history dating back to decades since the 
origin of statistical mechanics. It includes many different subjects and applications 
ranging from irregular motion of particles in homogeneous or disordered media and 
osmosis [H E], to standard and anomalous transport [31 H] of heat, mass and charges in 
materials, from coalescence of passive scalars [5] in highly turbulent fluxes to biological 
and ecological investigations on animal dispersal [ElUJ. 

The basis of the theory of diffusion is the random walk (RW) or random flight 
where a particle undergoes successive random displacements from its initial positions 
performing an irregular motion. The resulting character of the particle trajectory 
emerges from the statistical properties of these displacements. The wide applicability 
of the RW to natural phenomena relies just on the possibility to introduce appropriate 
generalizations on probabilistic nature of the displacements. One of the straightforward 
generalization is realized by introducing correlations in the displacements so to obtain 
the so called correlated random walks (CRW) [8]. This possibility extends also to 
a set of particles distributed in space leading to the definition of spatially correlated 
random walks. In this case one can wonder what kind of particle distribution emerges 
from reiterated displacements of the particles and how its properties can be directly 
inferred from the knowledge of the statistical correlation of the displacements P, [TO] . 
This approach is usually adopted to model the motion of a system of particles in a 
disordered environment where the randomness change in time. The study of this kind 
of models is strictly related to the mathematical theory of the so-called point-processes 
[TT] . i.e., stochastic point-wise particle distributions. The theory which investigates 
the geometrical and statistical properties of the stochastic processes generating particle 
distributions constitutes a very important field of research with important physical and 
cross-disciplinary applications [12] . 

In this paper we derive different results on the physics of diffusion from the exact 
transformation [13] of the two-point correlations of a particle distribution subjected to 
a single step stochastic deformation. More precisely, we deal with the transformation 
that a stochastic displacement field induces on given particle distribution once it is 
applied to the particles. Statistically independence between the particle distribution and 
the displacement field is assumed, but arbitrary displacement-displacement correlations 
and initial density- density correlations can be present. The main purpose of this study 
was originally to understand, through the exact formula presented for one a two-point 
particle correlations in [13], the discretization and finite size effects in the preparation 
of the initial conditions in cosmological ra— body [HI |15] simulations which are usually 
built by applying a suitable stochastic deformation field to a regular lattice or other 
very uniform configurations of particles with equal mass. In the present context, 
instead, we use the same approach to derive by reiteration all the famous equations 
of ordinary diffusion, super and sub-diffusion. We finally consider the most complex 
case of clustering and coalescence of particles observed in fully developed turbulence 
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due to spatially correlated diffusion induced by turbulent ffows on passive pollutants. 

2. From one step displacement to spatially correlated random walks 

In this section we recall first some notations for the one and two point correlation 
function of a spatially homogeneous point process, then we briefly summarize the main 
results of Ref. [T3] . 

Let us consider a d— dimensional spatial distribution of particles with equal 
unitary mass (i.e. a so called "point process") in a volume V which we assume to 
coincide asymptotically with M'^. The microscopic density is by definition 

N 

= ^5(X-Xi) , 

i=l 

where the limit is taken keeping fixed N/V = no, and Xj indicates the position of the i^^ 
particle. The system is supposed to be statistically homogeneous and uniform on large 
scales, therefore the average density is well defined and positive: (n(x)) = no > 0, where 
the average (..) runs over the ensemble of realizations of the point process. Moreover, 
we define the normalized connected two-point correlation function, as 

(n(xo)n(xo + x)) -^W , /,/ n /1^ 

The covariance function h{x.) is the off-diagonal two-point correlation function and gives 
the true correlation between different particles. The power spectrum of the particle 
distribution is defined as the limit 

S(k) = Ito - 2.„„^(k) , 

where ?7,(k; V) = jy d^x exp(— ik ■ x)?t,(x) is the Fourier transform (FT) of the density 
on the finite volume V . Because of the hypothesis on statistical homogeneity, 5'(k) is 
simply given by the Fourier transform (FT) of ^(x) multiplied by no: 

S(k)=no^[e(x)] = l+noMk), (2) 

with /i(k) = jF[/i(x)] and = / d'^'x e~''^'^{..) being the usual d— dimensional infinite 
volume Fourier transform. It is clear that for a statistically homogeneous particle 
distribution, S'(k) and ^(x) contains the same information. 

When a statistically homogeneous and arbitrarily correlated stochastic displace- 
ment field u(x) is applied to the particle distribution, each particle moves from its old 
position Xj to the new one Xj + u(xj). We assume that the displacement field and the 
particle positions are statistically independent. We are interested in the change of the 
two-point correlation properties of the particle distribution under the effect of the dis- 
placements. The complete statistics of u(x) is given by a probability density functional 
P[u(x)] giving the statistical weight of each realization of the stochastic field. However, 
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in our hypotheses, the behavior of the two-point correlation function or the power spec- 
trum under the system displacements [13] is only determined through the knowledge of 
the probability density function (PDF) 

0(w;x)= / [ d'^ud'^vf{u,v;x)5{w-u + v). (3) 



that two particles, separated by the vector distance x, undergo a relative displacement 
w. In Eq. ([3]) we used the properties that 0(w; x) is, in turn, related to the joint PDF 
/(u, v; x) that two particles separated by the vector x perform the displacements u 
and V respectively. If we denote by ^^(k) and .^m(x) = [6{x.)/no + hmi^x.)] respectively 
the power spectrum and the two-point correlation function of the particle distribution 
before the application of the displacements and by S'/(k) and ^/(x) = [5(x)/no + /i/(x)] 
the corresponding quantities after the displacements, we can write [13] the equation: 



Sfik) = l-J A0(k,q)+no J d'^a;e-'^''0(k;x)[l+e,„(x)]-(27r)'^no5(k),(4) 



where 



0(k;x) = j d'^we-'^-'^(l){w;x) 

is the characteristic functioiJil of the random displacement w, and 0(k, q) = 
J d'^x e~'^'^'^(f){\i.;x.). The only hypotheses for the validity of Eq. (jlj) are: (i) spatial 
homogeneity of both particle distribution and displacement field, (ii) statistical 
independence between the particle positions and the displacement field. Notice that 
no "small displacements" approximation or special properties of the displacement 
correlations are required. Let us define 



G^^(x) = m^(xo)m^(xo + x) /i, z/ = l,...,rf (5) 

the displacement-displacement correlation function which is a symmetric tensor of rank 
2, whose FT, is a non-negative definite symmetric tensor for all the fc-vectors. With (..) 
indicating the average over 'P[u(x)]. For symmetric distributions 'P[u(x)] = P[— u(x)], 
the positive value G^^{0) represents the variance of the fi-th component of u(x) at 
any point. Moreover in the case of a Gaussian field, the function G'^i,(x) determines 
completely the probabilistic properties of the field [T2] . 

The purpose of this paper is to deduce from Eq. (jlj) the time-evolution of the 
two-point correlation function and power spectrum of a particle distribution in the 
limit that each particle performs a Brownian trajectory, but the motions of different 
particles can be arbitrarily spatially correlated. In other words, we consider an assigned 
displacement field statistics, /(u, v;x), at every time step At, but no time correlation 
between consecutive time-steps. Then we take At — » in such a way to have a well 
defined diffusional limit. Two cases have to be basically distinguished: 

I From Eq. ([3]), it is immediate to verify also that (/)(k; x) ~ /(k, — k; x) and 0(k, q) ~ /(k, — k; q) where 
respectively /(k,k';x) = // /(u, v; x)e-*('' "+*''-^) and /(k,-k;q) = / d-^a; /(k, -k; x)e-*'i ''. 
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(i) The field u(x) is a spatially uncorrelated stochastic process at each time step, i.e., 
u(x) and u(y) are completely independent if x 7^ y. In this case. 



/(u, v;x) 



p{u)p{\) if X 7^ 
(5(u — v)p(u) if X = 



where p(u) indicates the single displacement PDF flB]. Basically, the displacement 
field is simply a white noise both in space and tim^. Moreover /(u, v; x) is a 
discontinuous function at x = 0, that is G^i^(x) = for x 7^ 0, whereas G^,^{0) 7^ 
and in particular G^u{0) = S^uU^/d > when assuming that p{u) = p{u). Then, 
we can simply show that 0(k; x) = |p(k)p where p(k) = jF[p(u)], and, because of 
the discontinuity in x = 0, J d'^q(f)(k; q) = |p(k)p 7^ 0(k; 0) = 1. This implies that 
Eq. dl]) reads 

sKk) = i+ip(k)n5..(k)-i]. (6) 

Note that this relation is local in k, i.e., each k— mode of the particle density evolves 
independently one of each other. 

(ii) u(x) is a real correlated and continuous stochastic process. In this case [16], it is 
well known that G^^{x) is a continuous function of x and in the limit when x — »• 0, 
it approaches with continuity its value G^u{0) [16]. In other words, /(u, v;x) is 
continuous in x and 

lim /(u, v; x) = 6{u — v)p(u) . 

X— >0 

this implies that J d'^q (pCk; q) = 0(k; 0) = 1 [13], and therefore Eq. (jlj) becomes 

Sf{k) = noj d'^xe-^'^-^^ik; x)[l + ^^.(x)] - (27r)'^no5(k) . (7) 

Unlike Eq. (I6l), this equation is nonlocal in k due to the presence of displacement- 
displacement spatial correlations which couples different modes of the particle 
density before and after the application of the displacement field. 

We can now study the spatial diffusion of a particle distribution in which at 
each time steps the particles move under one of the stochastic displacement field just 
described above. We will assume for simplicity, that the time can be discretized in time- 
steps of size At at which the particle distribution is displaced by a realization of the 
field u(x, t); moreover different time-steps are supposed to be statistically independent. 



3. Spatially uncorrelated displacements 

In order to illustrate the general formalism to derive a continuous time equation from 
the discrete one (single step), it is instructive to consider the simplest case of spatially 

§ Apart from the possible correlations in d > 1 between the different components of the displacement 
u(x) at a single point x. However, if we consider the case that p(u) depends only on u = |u| (i.e., the 
displacement field is isotropic), perpendicular displacements are uncorrelated. 
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and temporally uncorrelated displacement field. As we show below, it corresponds to 
the homogeneous diffusion equations, either standard or fractional depending whether 
the variance of displacements is finite or infinite. This discussion is useful in view of the 
more interesting case of random walks generated by spatially correlated displacements. 
As clarified above, all the statistics of the field is contained in the one-displacement 
PDF p{u). Let us firstly consider the statistically isotropic case for u: p{u) = p{u), 
implying that p(k) = p{k). We have to distinguish the two cases of finite and infinite 
variance u"^. In the former, each particle performs an ordinary (i— dimensional random 
walk and the paths of different particles are independent one of each other, while in the 
latter, each particle undergoes a (i— dimensional Levy walk independently of the others. 



3.1. The continuous time limit 

For statistical isotropic displacements, we can expand at small k as 

p{k) = 1-Bk" + o{k'') , (8) 

where a = 2 and B = /2d when v? is finite, while < a < 2 when v? diverges and 
p{u) ~ Au~'^°'^''^ at large m, with B proportional to the amplitude of the tails A > 0. In 
the first case each particle undergoes an independent standard random walk while, in 
the second, an independent Levy flight [17J. A well deflned diffusional continuous time 
limit is attained by requiring that 2B = DAt where D > is a constant independent 
of the time step At. Therefore in the limit At ^ 0, the substitution of Eq. ([8]) into 
Eq. §i) leads, for A;" < 2 /{DAt), to 

dtS{k,t) = Dk''[l-S{k,t)] , (9) 

whose solution is 

S{k, t) = [SM - 1] exp(-DA;"t) + 1 . (10) 

The approach of ^(k, t) to homogeneous Poisson power spectrum S'(k) = 1 is 
exponentially fast, and each k— mode of the two-point correlation relaxes with a rate 
Dk" (0 < a < 2). We see immediately that the smaller a the faster the approach to the 
completely uncorrelated stationary state. In other words inhomogeneities diffuse (a = 2) 
or super-diffuse (0 < a < 2) until reaching the uniform stationary state rs(x) = no. 
Equation ([9]) can be recast in a more familiar form by rewriting Eq. (|TOi) for h(k,t): 

dth{k,t) = -Dk''h{k,t) . (11) 

and then taking the inverse FT. For simplicity, let us consider a one dimensional system. 
By taking the FT of Eq. f|TT]) we can write 

dth{x,t) = Dd;^h{x,t) , (12) 

where is for a = 2 the usual Laplacian in d = 1, giving the celebrated standard 
diffusion equation, while for < a < 2 is the fractional derivative [H] of order a 
in X. Equation (fT2l) is for < a < 2 the well known equation of spatially fractional 
diffusion [181 IlH] and describes how the connected two-point correlations super-diffuse 
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in space toward an uncorrelated Poisson steady state. The equation (1121) takes a more 
familiar form when written for the average conditional density r{x,t): 

dtT{x,t) = Dd^r{x,t) , (13) 

which is the fractional diffusion equation for the density of particles seen in average by 
a generic particle of the system. It belongs to a larger family of fractional diffusion 
equations whose general form is 

d^r{x,t) = Dd;^r{x,t) (14) 

with < (3 < 1 and < a < 2. This kind of equation is encountered in the description of 
key aspects of anomalous transport as for instance that occurring on disordered peculiar 
structures (fractal supports) or due to spatial non locality [HI [20]. More precisely the 
fractional diffusion equation ([T3j) can be obtained in the context of continuous time 
random walks (CTRW) [21] under the assumption that the joint PDF tp{u,t) to make 
a step of size u in the time interval [t,t + dt) factorizes as 0(t)p(u). In order to obtain 
/? = 1 it is necessary that the mean value of 0(t) is finite, and to obtain « = 2 we need 
a finite variance for p{u). Otherwise we have respectively (3 < 1 and a < 2 depending 
on the power law tails of the two functions 0(t) and p{u). Indeed we have obtained 
Eq. f|T3|) via the choice (j){t) = 6{t — At) (with At 0^) which has finite mean value 
equal to At itself. 



4. Spatially correlated displacements 

We now turn to the case where the stochastic displacements, acting on the particles 
at each time-step, are generated through the realizations of a continuous and spatially 
stationary correlated stochastic field u(x) defined by a probability density functional 
7^[u(x)]. For instance, in d = 1 we can consider a Gaussian stochastic field defined by 
the probability density functional 



V[u{x)] ~ exp 




dx dyu{x)K{x — y)u{y) 



with K{s) being the positive definite and continuous correlation kernel. Since we assume 
again no time correlation between successive realizations of the field u(x), at each 
time-step the power spectrum of the particle distribution evolves according to Eq. (^^. 
In principle, a time dependent PDF 'P[u(x);t] can be also considered, referring to a 
displacement field whose correlation properties depend on time. Here we limit the 
discussion to time independent functionals V, focusing on the case of a dimensional 
Gaussian displacement field with finite variance u"^ = X])!=i^mm(0) < +oo, where 
the displacement-displacement correlation matrix G^,^(x) has been defined in Sec. |2j 
However, as shown explicitly below, the validity of the evolution equations we derive 
is not restricted to this case. If the displacement field is Gaussian its characteristic 
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0(k;x)=exp|-^A;^A;4G^,(O)-G'^,(x)]| , (15) 

where A;^ indicates the /i*^ component of k. 

In analogy with the simpler case of uncorrelated displacements, we have to assume 

G'^^(x) = Atc^^{x) 

for the existence of a smooth time limit, with c^;y(x) independent on At in order to obtain 
the correct diffusional processes in the limit At — > 0. In fact the above condition implies 
that the variance of each component /i of the field satisfies G^^{0) = |u^(x, t)^ ~ At at 
any time t. Therefore, if u(x, nAt) is the displacement in the point x at time t = nAt, 
we can write: 



u(x,nAt) = (16) 
M^(x, nAt)M^(x', n'At) = Atc^^(x - x')(5„„/ , (17) 

thus, in the limit At ^ the field ?7(x, t) = u(x, t) / y/At becomes a spatially correlated 



and temporally delta-correlated noise for the motion of particles. In a hydrodynamic 
analogy, in which the displacement field is interpreted as a turbulent main flow advecting 
the particles of a passive pollutant, the quantities d^i,{'x) = [c^,y(0) — c^,y(x)] are usually 
called [22] the structure functions of the flow. 

In the limit of small At, we can expand Eq. ( fTSjl 

l,d 

0(k; x) = 1 - At ^ k^Kd^A^) + o(At) . (18) 
When plugged into Eq. ([7]) this gives 

S{k, t + At) = 5(k, t) + At ^ k^K [noc^uO^) - c^.(0)5(k, t)] + (19) 



+ Atf^k^k,J A^S(k-q,t)c^,(q) + o(At), 



where c^i,(k) = J^[c^^{:>i)] is the re-normalized power spectrum tensor of the 
instantaneous displacement field. The limit At — *• on Equation f[T^ yields 

^''^ r f d'^q - 

(k)-c,.(0)5(k,t) + j ^5(k-q,t) 



which in terms of FT of the covariance function takes the same forms 

dMKt) = J2k,k (k) - c^^(0)/i(k,t) + / — -y/i(k- q,t)c^.^(q) 



■(2i; 
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l,d 

dth{x, t) = J2 ^lu MM-(x)/i(x, t) - c^,(x)] . (22) 

IJ,,U 

Consequently, for r(x) = no[l + /i(x)] we obtain 

l,d 

d,T{^,t) = J2dlAd,u{^n^,t)] . (23) 

It is noteworthy that if P(x, t|xo,to) is the PDF of the separation x at time t between 
an arbitrary pair of particles of the distribution, given their initial distance xq at time 
to, we can write 



r(x, ^) = j d'^^o ^(x, i|xo, to)r(xo, to) . 



It is straightforward to verify that the transition probability P(x, t|xo,to) also satisfies 
Eq. and for this reason P(x, t|xo, to) is called the propagator of the diffusion operator 
defined by Eq. fl23|) . On the other hand Eq. fl23|) is the Fokker-Planck (FP) equation 
associated to the stochastic Langevin equation (LE) in Ito representation [23j 

x(t) = w(t) (24) 

describing the time evolution of the two particle separation x(t), where w(t) is a 
Gaussian noise with the following one and two-time correlation properties: 



w^{t) = 

/i,z/ = l,...,(i (25) 

w^{t)w,{t') = 2c/^,(x)(5(t - t') 

Therefore Eqs. f l23|) and fl24|) are equivalent. 

In the hyper-isotropic condition, corresponding to the choice (i/^i/(x) = 5^ud{x) [i.e., 
c^i,(x) = 5^yc{x)], Eq. (1251) can be rewritten as 

9ir(x,t) = V2[ci(a;)r(x,t)]. (26) 

Note that all these results are not restricted to the case of Gaussian displacement fields. 
In fact, as shown in [13], if G^yi^x) = u^{xo)uu{xo + x) is finite for a; = 0, the small k 
expansion 

l,d 

0(k; x) = 1 - ^ k,K[G,,iO) - G^,(x)] + o{k^) 
is always valid. 

It is easy to prove, but important to note, that when the displacement field 
u(x, t) can be decomposed into two independent components: a spatially correlated 
field Ui (x, t) characterized by Eq. (ITTj) and a spatially uncorrelated and statistically 
isotropic field U2(x, t) of variance u| = d x DAt, Eq. (!23l) takes the form: 

l,d 

dtT{x,t) = DV2r(x,t) + ^2^1^ [rf^.(x)r(x,t)] . (27) 
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That is, the motion of the set of particles results in the superposition of a standard 
diffusion, first term in Eq. (j27|) . with a spatially correlated diffusion, second term. 



Equations of the type fl2711 are the generalized diffusion equations which are usually 
encountered in the context of turbulent transport of passive scalars (i.e., pollutant) 
[5], where ?t,(x, t) is the density of the passive particles advected by the velocity field 
v(x, t) = u(x, t)/At of the synthetic turbulent flow. It is important to note that for 
x — > 0, the term DV^r(x, t) of Eq. fl27|) dominates because dfj,^{0) = 0, therefore the 
small scale motion occurs via standard diffusion. The interesting case is obtained when 
D — > 0, and consequently x = becomes a singularity of Eq. (!27|l |23j . In particular in 
turbulence a complete solution has been given for the scale free case in which 

df,^{x.) = ax^5^^ + hx^~'^x^Xy , (28) 

with < ,^ < 2, a and h constants such that the tensor c^!y(x) has a positive definite FT. 
The one dimensional case is recovered by putting d = 1 and 6 = 0. The class of models 
defined by the structure tensor (!28l) is referred to as generalized Kraichnan ensemble 
[211 [5] and a complete classification of their solutions, in terms of the singular behavior 
around x = 0, has been given in Ref. [25l[26]. For < < 2, it consists in three possible 
behaviors: 

(i) For {h/a) > [d — 2) and at the same time (^ — l){b/a) > [d — C,), the only 
possible solution of Eq. fl23l) is such that different particles coalesce in finite time 
and no stationary state exists. In practice pair of particles collide in finite time 
with vanishing relative velocity and therefore remains attached for the rest of 
the dynamics. In this case r(x) ~ a{t)6{x.) + (3{t)x'^~^~'^ with a time increasing 
coefficient a;(t), which signals the coalescence phenomenon; 

(ii) For (^ — l){h/a) < {d — ^) and at the same time {hja) < {d — 2) the only possible 
solution is such that the probability of finding more than one particle in a single 
spatial point is zero at all time. In this phase particles form only clusters and the 
diffusion of particles converges to a stationary state in which r(x) ~ x""^ where 
7 = ^ + 6(d - l)/(a + 6); 

(iii) Finally for (b/a) > {d — 2) and (^ — l)(6/a) < {d — C,) particles can collide at finite 
time, but with non-zero relative velocity. This obliges to fix a boundary condition 
at X = for the diffusion equation. For an absorbing boundary condition one has an 
effective behavior as in case (i) above. Instead for a reflecting boundary condition 
the effective behavior is similar to case (ii) above. Adopting a mixed boundary 
condition a composition of the two above behaviors appears (called "sticky" phase 
in [26]). 



The method used to get this classification is quite complex and consists in the theory of 
boundary conditions of elliptic operators. In the rest of the paper we do not enter the 
details of this rigorous analysis, but limit our study to the one dimensional case for a 
generic choice of the displacement-displacement correlation function through computer 
simulations and simple theoretical arguments. 



Diffusion, super- diffusion and coalescence from single step 



11 



4.1. The one dimensional spatially correlated diffusion 



In d = 1 Eq. ( !23l) becomes 



dtT{x,t) =dl[d{x)T{x,t)] 



(29) 



where d{x) = [c(0) — c(x)] is the structure function of the stochastic velocity field (i.e., 
displacements) with two-point correlation function c(x). It is important to note that the 
general properties of any correlation function c{x) constraints d{x) to have the small 
X behavior d{x) ^ ax^, with a > and < < 2. In this paper we do not treat 
the "smooth" case ^ = 2 as it has been solved rigorously elsewhere [10] and focus our 
analysis on the "rough" stochastic velocity fields for < ^ < 2. The aforementioned 
classification of the solutions of the Kraichnan ensemble immediately implies that for 
5 = and d = 1 in Eq. fl28l) . only the first and the third case are possible in one 
dimension around the singularity x = 0. In particular we have the first behavior for 
1 < ^ < 2, and the third for < ^ < 1. For ^ = 2 (when the velocity field is smooth at 
small scales) one can see that particles coalesce but in an infinite mean time p!0]. 

First of all we note that, if it exists, the only possible stationary solution of Eq. (!29ll 

is 



Ts{x), due to its definition of average conditional density, must converge to hq for 
|x| — * 00 where density- density correlations has to disappear. In order to decide around 
the acceptability or not of this stationary solution, it is necessary to study its small scale 
behavior. We know that by definition Ts{x) is acceptable only if it is integrable at small 
X. We see immediately that in = 1 it happens only for ^ < 1. Thus for ^ > 1, Eq. (!29ll 
admits no stationary solution in agreement with the aforementioned classification of 
solutions of the generalized Kraichnan's ensemble In fact Eq. fl30p corresponds 

to the stationary state which the family of solutions with simple clustering and no 
coalescence converge to, obtained by imposing a reflecting boundary condition at the 
singularity x = 0. This solution is analogous to the Poisson stationary correlation 
function Ts{x) = uq for the ordinary diffusion equation to which corresponds the well 
known propagator 



satisfying the reflecting condition at x = and which shows how two particles at initial 
vanishing distance spread when a reflecting condition is imposed about their collisions. 
An analogous propagator can be found also in relation to the stationary state (|5UI) of 
our more complex diffusion Eq. fl29|) . It can be found by looking for a scaling solution 
of the form P{x, t\0, t) = t~^ f{x/t^). Plugging this scaling form into Eq. ( !29l) . one finds 



r^(x) = no 



c(0) 



(30) 





(31) 




(32) 
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where C > is the normahzation constant. Also Eq. (1321) describes the spreading of 
a pair of particles at initial vanishing distance with reflecting condition at x = 0. In 
fact dx[x^P{x,t\0,t)]\x=o = 0. Moreover it becomes the ordinary diffusing Gaussian 
for ^ — i> 0. Let us now analyze in detail Eq. (!30|) . First of all we note that 
for all X the covariance function h{x) = [T{x) — uol/no has the same sign of the 
displacement correlation function c{x), i.e., those scales at which displacements are 
positively (negatively) correlated asymptotically become also positively (negatively) 
correlated scales for the density of particles. Moreover Ts{x) for large x approaches 
the average density Uq in the following way 

c(x) 



c(0) 



(^xj ^ no 1 + 
i.e., h{x) is 

c{x) 

" W) • 

In other words the iterated displacement field injects exactly its large scale correlations 
in the particle system. This is interesting because, as one can check by expanding 
Eq. (171) at small k [I3], the displacement field injects, in a single step, only a large scale 
contribution to the power spectrum S{k) [or h{k)] of order k'^c{k)At which vanishes k"^ 
times faster than c{k). Finally, at small x we have Ts{x) ~ x~^, meaning that particles 
forms, at small scales, clusters with fractal dimension D = {1 — 

We now turn to the problem of what happens for > 1 and how the other "non- 
reflecting" solutions for ^ < 1 behave. This is a more difficult task as the singularity 
at X = of Eq. (l29l) generates a coalescence dynamics which can not be described 
through only smooth functions. In order to study this case, we adopt a sort of mean 
field approximation in an appropriately transformed LE (in Ito representation) for the 
separation x between a pair of particles associated to the FP Eq. fl29l) . We suppose that 
at initial time t = we have xq = x{0) > 0. At sufficiently small x Eq. ( 12^ . for our one 
dimensional case, can be written as 

x{t) = A[x{t)f'r^it) , (33) 



where A = \'2a > and ri(t) is a white Gaussian noise such that ri{t) = and 
ri{t)r]{t') = 6(t — t'). Let us now apply the following change of variables to Eq. ( l33i) : 



I.e., 

2/(2-5) 

(35) 



X 



By using the rule of change of variables in the Ito representation [23] , we get the LE for 

yit) 

yit) 
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where C = ^/(4 — 2^). If we suppose that the initial value y{0) > is sufficiently small, 
the noise ri{t) can be neglected obtaining 

C 



yit) 



whose solution is 



yit) = ^ym-^^t, (36) 

l-{/2 

where y{0) = from Eq. ([M]). This equation implies that at t* = [2 - Oy^(O)/^ 

we get y{t*) = and therefore x{t*) = 0, i.e. the two particles collide. Moreover the 
velocity x{t), in the approximation where we neglect the noise, is 

(C-l)/{2-0 



. , , dx dy 

x(t\ = ■ — 

dy dt 



Therefore x{t*) = for 1 < ,^ < 2, while x{t*) —>■ — oo for ,^ < 1. Consequently, for 
1 < ^ < 2, independently of the boundary condition at a; = 0, once x{t) vanishes 
it keeps this value forever, i.e., for 1 < ^ < 2 this coalescing solution is the only 
possible solution as predicted by the general classification given above for the Kraichnan 
ensemble. Instead for ^ < 1, since the velocity x{t) is non-zero, when x{t) = 0, it is 
necessary to choose by hand the boundary condition to fix the kind of solution. For 
reflecting boundary condition we have the above presented solution converging to the 
stationary state Ts{x) given in Eq. ( l30l) . while for absorbing boundary condition no 
stationary solution is reached and we have a behavior similar to that for ^ > 1. We now 
derive the small x behavior of the propagator P{x, t\xQ, 0) with sufficiently small Xq > 
from Eq. (l36il valid for both cases < ^ < 1 with absorbing boundary condition and for 



1 < ^ < 2. Equation fl36l) says that at time t all pair of particles with ?/^(0) < 2Z|t, i.e., 

with initial separation xq < = Xmaxit) have already collided at time t. 

For 1 < ^ < 2 such pairs of particles always coalesce, as x{t) = when x{t) = 0. Instead 
for < ^ < 1 they coalesce only if we impose the absorbing condition (i.e. completely 
inelastic collision) at a; = 0. In both cases we can say that P{x, t\xo, 0) has developed a 
singular contribution m{t)6{x) at x = where 

m{t) = / (ixop(xo) , 



where p{xo) is the PDF of the initial pair distance Xq. If Ts{x, 0) = uq (i.e. initial particle 
distribution is a homogeneous random Poisson one) for small Xq (i.e. Xq <^ I/'t-o) we 
have p{xo) ~ no, from which we obtain m(t) ~ noXmax{t) ~ t^^^'^~^\ On the other 
side for all pairs with xq > Xmaxit) we can develop the following argument. Since in 
the present approximation the dynamics is deterministic one can derive P{x,t\xo,0) 
at a; > directly by a simple change of variable in p{xo). More precisely let us call 
Xt = x(t) and P{x, t\xo, 0) = Pt{xt); we can write by conservation of probability: 

Ptixt) = p{xo)^ (37) 
dxt 



Diffusion, super- diffusion and coalescence from single step 



14 



and the relation between Xt and Xq is given by Eqs. and (15^ . By considering 

again ^(xo) — ?t-o, it is simple to derive for small x > that 

Note that lim^j-^Q pt{x) = for ^ < 1 as it has to be for an absorbing condition at x = 0. 
Finally at small x and sufficiently small t, if r(x — ^> 0, 0) = hq, the expression 

T{x, t) ~ P{x, t\xo, t) = + ct«-^)/(2-«)xi-« , (38) 

represents the solution for 1 < ^ < 2, while it is a solution for the case < ^ < 1, only 
when considering absorbing conditions aX x = 0" . With h and c positive constants. This 
completely agrees with what found in [25l [26] . 



4-.2. Numerical results 

In the light of the previous one dimensional analysis, we discuss an important example 
m. d = 1 where u{x) is a spatially correlated Gaussian displacement field with short 
range correlation function 



c[x) 









L 


X 






Xq 





< ^ < 2 . (39) 



where Xq determines the correlation length. We recall that according to the Wiener- 
Khinchin theorem [H], a given c(a;), such that c(0) > 0, is a well defined correlation 
function of a continuous stochastic field if and only if c{x) is continuous for all x, c{k) > 
for any k, and finally dk c{k) < +oo. All these conditions are satisfied by Eq. (!39|) . 
Varying ^ between and 2 allows us to study numerically all the phases above described 
theoretically: 

(i) for < ^ < 1, in case of refiecting boundary conditions on particle-particle collisions, 
the particle diffusion admits the stationary solution Eq. (l30l) implying a small x behavior 



[X) ^ no 



(40) 



X 
Xq 

which describes a fractal distribution with dimension Df = 1 — ^ for x < Xq, and exhibits 
a cross over to a uniform particle distribution with average density uq for x > xq. For 

X ^ Xq instead correlations decay exponentially fast as hs{x) ~ exp 

(ii) for 1 < ^ < 2, the system admits no stationary behavior and the unique solution 
corresponds to particle coalescence as described by Eq. (155]) . 

We run computer simulations to check these behaviors of particle distributions 
under the effects of a repeated application of the displacement field with correlation 
( l39l) . The numerical implementation of the dynamics requires the generation of an 
array of correlated Gaussian random variables {ui} with the prescribed correlator 
(uiUj) ~ c(xj — Xj) where c(x) is given by Eq. (|39l) and Xi is the position of the i*^ 
particle {i = 1, ...,N). We used two methods: 
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time 



Figure 1. Detail of the space-time pattern generated by the trajectories of = 512 
random walkers undergoing the evolution Eq. (|41|) with time step h = 10^^ and 
displacements that are not time-correlated but have a spatial correlation (j39|l . with 
xq = 3. Left panel corresponds to ^ = 0.5 showing simple clustering and right refers 
to ^ = 1.5, for which a coalescence regime occurs. 



(i) The first one considers the Cholesky decomposition method [27] for the covariance 
matrix. According to this algorithm, an array {ui, ....,un} of N correlated Gaussian 
variables with correlation matrix C is obtained from the set of {wi,....,wn} 
independent random Gaussian variables with zero mean and unitary variance by 
applying the linear transformation u = Aw, where A is the lower diagonal matrix 
A such that C = AA* (Cholesky decomposition of the matrix C). 

(ii) The second method makes use of discrete Fast Fourier Transform (FFT) to generate 
the stationary Gaussian field {u{xn)} over a grid of sites x„ with n = 1, ...,J\f. The 
values of the field on the site Xn is expressed as the Fourier sum 

k 

with k = 27rm/A/'and m = 0, 1, to avoid aliasing. Note that z{k) = z*{—k) 
to ensure that u{xn) is real. The choice z{k) = [a{k) + iP{k)]/2 with a{k) and P{k) 
independent real Gaussian variables of zero mean and variance Afc{k)/2, guarantees 
the set of random variables representing the discretized version of the field 

to have the correct power spectrum c{k). 

Once the fields is generated the position of a particle i is updated according to the Euler 
scheme 

Xi{t + h) = Xi{t) + Vhu{xi) (41) 

with time step h. When the Gaussian displacement field is generated through FFT, 
there is clearly a problem associated to discretization due to the grid of step Ax where 
the FFT is computed. We assign a particle the displacement m„ if its position at time t 
falls in the n-th bin [nAx, {n-\-l)Ax] determined by the grid. The choice of the time step 
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is such that h{u{xiY) — i'^ order to sample the maximal resolution scale allowed 
by the discretization. 

Spatio-temporal patterns obtained via simulations of a system of particles 
subjected to displacements of correlation c{x) with xq = 3 are shown in Fig. [T]for cases 
^ = 0.5 and ^ = 1.5. The particle are initially distributed uniformely in a 1-dimensional 
simulation box with density p = 1, periodic boundary conditions are applied at the ends 
of the box. For ^ = 0.5 (left panel) and assuming elastic collisions (a reflecting boundary 
condition when particle trajectories intersect one another), the system exhibits simple 
particle clustering toward a stationary state. Instead for ^ = 1.5 (right panel), particle 
trajectories coalesce more and more in time and no stationary regime is actually reached. 
We measured during each run the density correlation function r(x), as the histogram of 
the relative particle distance \xi — Xj \ and the results, averaged over several independent 
runs starting from the uniform particle distributions, are plotted in Figs. [2] and [3l for 
^ = 0.5,0.95, 1.5 respectively. The stationary state for the case ^ = 0.5 coincides with 
that described by theoretical Eq. (l30ll which is represented in Fig. [2] (together with the 
one of the case ^ = 0.95 which is in the same class of behavior). The discrepancy in the 
amplitudes between numerical and theoretical stationary Ts{x) is due to the following 
finite size effect: in simulations of of particles in a volume V , the average conditional 
density T{x) is subjected to the integral constraint V{x)dx = {N — 1). Consequently, 
if the dynamics develops a strong positive density correlation at small scale, an artificial 
negative correlation must appear at larger scales. The typical behavior of r{x,t) for 
^ = 1.5 is represented in Fig. [3l It shows correctly the small x behavior proportional to 
x^-i predicted by theoretical results. 

5. Conclusions and Discussion 

The evolution of several many particle system can be described and also generated via 
the iterated application on the particles of suitable displacement field which determines 
their trajectories. In this context, the question we addressed concerns the connection 
between the statistical properties of the displacement field at a single time and the 
spatial correlations that arises in the particle distribution during the evolution. We 
have shown how to derive the partial differential equations describing the continuous 
time evolution of the two-point correlation function of the particle density under the 
iterative application of a stochastic displacement field with no temporal memory from 
the corresponding exact single step evolution equation. This continuous time equations 
are of a diffusion type and describe simple, fractional, or spatially correlated diffusion 
of density fluctuations, depending only on the two-point statistical properties of the 
elementary displacement field. Simple and fractional diffusion occurs in absence of 
spatial correlations of displacements. Which of the two regime arises depends only on the 
finiteness of the elementary displacement variance. Spatially correlated displacements, 
instead, determine a Fokker-Planck (FP) equation for density correlations equivalent to 
a simple Langevin equation whose multiplicative noise is defined again by the two-point 
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Figure 2. Log-log behavior of tlie average conditional particle density r(a;) as function 
of the separation x generated by the iterations of Eq. (|4T|) with time step h = 10^'* 
and a displacement field with correlation (j39|) defined by parameters xq = 3, = 0.5 
(Left) and ^ — 0.95 (Right). T{x) is computed as the histogram of the intcr-particle 
distance binned exponentially in 128 intervals. The number of particles is = 256 and 
simulation data are the results of the average over 1500 independent runs. The small 
scale decay is expected to be a power law with exponent ^. The solid line indicates the 
result of a power law fitting with exponents —0.5, —0.95 respectively and the dashed 
lines refer to Eq. (|40|) . 



r 




X 



Figure 3. Log-log plot of the average conditional density r(a;, t) as function of the 
separation x at time t = 8 , for ^ = 1.5, xq — 3 and = 64 particles, averaged over 
2200 independent runs. Integration time step in Eq. (j4T|) amounts to h = 10~^. The 
small scale decay is fitted by a power law (solid line) with exponent 0.52 in consistent 
with the theoretical value (1 — ^) = 0.5. The inset shows the time increasing behavior 
of the coefficient of the 5 function contribution to Eq. (j38p due to particle coalescence. 
The solid line indicates the mean field scaling behavior ~ 
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displacement correlation function. Such a FP equation predicts a rich phenomenology, 
ranging from simple particle clustering to coalescence. These two regimes depends 
uniquely on the small scale behavior of displacement correlations. 

We characterized and classified these regimes in terms of the the solutions T{x,t) of 
the FP equation for the density correlations which depend on the small scale statistical 
properties of the displacements fields. This classification can be given in terms of the 
theory of boundary conditions of FP equation [28j. Our problem, indeed, presents 
in general a singularity at the origin corresponding to the vanishing of the structure 
function of the velocity field and which has to be treated as a possible additional 
boundary [23]. By developing a simple theoretical approach in = 1, we have compared 
the results with the solution classification for the Kraichnan ensemble for d— dimensional 
turbulence [25] finding a perfect agreement. 

Computer simulations, implementing the evolution of a one dimensional system of 
particles driven by a Gaussian correlated displacement field, confirm the presence of the 
clustering and coalescence regimes in agreement with our theoretical predictions. 

It is noteworthy that the problem we have considered is different from the motion of 
particles in random potentials for which in general different FP equations hold [28] . This 
point is simply clarified by the observation that, in our model, for any spatial correlation 
of the displacement field, any single particle performs a simple Brownian motion, and the 
effect of the two-point displacement correlations amounts to correlating the Brownian 
trajectories of different particles. On the contrary in presence of an external random 
potential even the motion of a single particle is far from being simply Brownian. 

It is also interesting to discuss the analogies with the random one dimensional 
coagulation process A + A A [29j, where particles perform mutually independent 
random walks and when they collide, one of them is removed according to an assigned 
rule. The particles are considered to be "sticky" as they undergo perfectly inelastic 
collision and coalesce in a single particle independently of the relative velocities. In our 
case on the contrary, for 1 < ^ < 2, particles are not "sticky" and their coalescence is only 
a consequence of the small scale correlation properties of the displacement field which 
imply particle-particle collisions occurring at vanishing relative velocity. Therefore, even 
though we would impose perfectly elastic collisions, particles would coalesce the same. 
The process A + A ^ A is recovered in the limit ^ — > (implying collisions at non-zero 
relative velocity with probability one) and imposing inelastic collisions at interparticle 
contact. 

In summary, via simple physical arguments, we have generalized the result of 
Ref. [To], concerning only the smooth case ^ = 2, to the case of < ^ < 2 where ^ 
is the exponent defining the behavior of displacement correlation c(x) at short scales: 
d{x) = [c(0) — c{x)] ~ x^. 
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